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ABSTRACT 

We present results of an extinction-CO line survey of the southeastern part of the California Molec¬ 
ular Cloud (CMC). Deep, wide-held, near-infrared images were used to construct a sensitive, relatively 
high resolution ( ~ 0.5 arc min) (NICEST) extinction map of the region. The same region was also 
surveyed in the ^^CO(2-l), ^^CO(2-l), C^®0(2-1) emission lines at the same angular resolution. These 
data were used to investigate the relation between the molecular gas, traced by CO emission lines, and 
the dust column density, traced by extinction, on spatial scales of 0.04 pc across the cloud. We found 
strong spatial variations in the abundances of ^^CO and C^®0 that were correlated with variations in 
gas temperature, consistent with temperature dependent CO depletion/desorption on dust grains. The 
^^CO to C^®0 abundance ratio was found to increase with decreasing extinction, suggesting selective 
photodissociation of C^®0 by the ambient UV radiation held. The effect is particularly pronounced 
in the vicinity of an embedded cluster where the UV radiation appears to have penetrated deeply 
(i.e.. Ay < 15 mag) into the cloud. We derived the cloud averaged X-factor to be <Xco> = 2.53 x 
10^° cm“^ (K km s“^)“^, a value somewhat higher than the Milky Way average. On sub-parsec scales 
we hnd there is no single empirical value of the ^^CO X-factor that can characterize the molecular 
gas in cold (Tk 15 K) cloud regions, with Xco oc Ay^^ for Ay > 3 magnitudes. However in regions 
containing relatively hot (Tex ^ 25 K) molecular gas we hnd a clear correlation between W(^^CO) 
and Ay over a large (3 < Ay < 25 mag) range of extinction. This results in a constant Xco = 1-5 x 
10^° cm“^ (K km s~ )“^ for the hot gas, a lower value than either the average for the CMC or the 
Milky Way. Overall we hnd an (inverse) correlation between Xco Tex in the cloud with Xco oc 
Tex~°'^. This correlation suggests that the global X-factor of a Giant Molecular Cloud (GMC) may 
depend on the relative amounts of hot gas contained within the cloud. 

Subject headings: stars: formation; ISM: clouds; (ISM:) dust, extinction; ISM: abundances 


1. INTRODUCTION 

The relationship between the dust and gas in a molecu¬ 
lar cloud is crucial in understanding the cloud properties 


and the star formation activities therein. Dust provides 
the most reliable a nd useful tracer of the total gas col¬ 


umn density (e.g., Bohlin et al. 1978 Rachford et al. 


2009 


Goodman et al. 20091, while molecules like ^'®CO 


and C^®0 are more useful as tracers of the local physi- 
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cal and chemical conditions (e.g., gas temperatures and 
velocities, CO column densities, etc.) A comparison be¬ 
tween them allows us to probe the variation of molecular 
gas abundances, i.e., [^^CO] and [C^®0] ([molecule] = 
N(molecule)/N(H) where N is column density), caused 
by such factors as photodissociation and chemical deple¬ 
tion/desorption in cold/warm environments. This helps 
to reveal the physical and chemical properties of the 
molecular cloud, which enables us to better understand 
the initial conditions of star formation and its poten¬ 
tial dependence on such factors as cloud chemistry and 
structure. Such variations have been noted in previous 


(Pineda et al. 20081 

, Taurus (Pineda et al. 

2010), and 

Orion ([Rippie et al. 

2013), made at relative! 

y high spa- 


of a combined CO-dust e xtinction study of the California 
molecular cloud (CMC, |Lada, Lombardi, & Alves|2009 1 
at even higher spatial resolution ( ~ U.U4 pc) in order to 
investigate in more detail the relation between CO gas 
and dust in this interesting nearby CMC. 

The CMC is a useful laboratory for such studies. In 
particular, it contains a massive filamentary structure in 
its southeast region that consists of both cold molecular 
gas with little star formation activity and hot molecular 
gas associated with an active star forming region con¬ 
taining the massive star LkHa 101 and its accompany¬ 
ing embedded cluster. This provides a physical environ¬ 
ment with large temperature and density variatio ns that 
are interesting to explore . Indeed, as reported in |Lada,| 


Lombardi, & Alves (2009), the CMC overall shows an or- 
der of magnitude lower star formation rate (SFR) than 
the Orion molecular cloud, even though these equally 
distant G MCs have simi lar mass, size, and morphology. 
Recently, Li et al. (2014) also found the ^®CO clumps in 
the CMC to be ch aracterized by similar kinematic s tates 
as those in Orion. Lada, Lombardi, & Alves (2009) pro¬ 
posed that the difference in Sk'Rs between the two clouds 
was a result of a difference in cloud structure, specifi¬ 
cally the amount of dense gas in each cloud. Moreover, 
this close connection between dense gas mass and the 
SFR appea rs to be a general physical property of galac¬ 
tic GMCs (Lada et al. 2010). It would be of interest to 
compare other properties of the CMC with similar ones 
in the Orion cloud to look for additional factors, such 
as cloud chemistry for example, that could contribute to 
the difference in the SFR between the two clouds. 

A combined CO-extinction study can also be used to 
investigate the CO X-factor (Xco = N(H 2 )/Wco(i-o)j 
where W is integrated intensity) which is widely used to 
derive molecular cloud masses from the CO(l-O) line flux, 


especially in external gala xies (see e.g. Sandstrom et al. 


2013 Bolatto et al. 2013). Despite being very optically 


thick, this CO rotational transition is still thought to be 
a good tracer for the total molecular gas mass, due to the 
effect of photon transport in gas with large velocity gradi¬ 
ents. However, variations of physical conditions can have 
an effect on the X-factor and, in many instances, render 
its application rather uncertain. Apparently, th ere is no 


universally valid X-factor (see recent review from Bolatto 
et al.|2013D . Consequently, investigating the dependence 


of the X-tactor on differing physical conditions is clearly 
of great interest in establishing practical guidelines for its 


application. Combined CO-extinction studies of nearby 
clouds can therefore provide potentially important in¬ 
sights into this issue. 

In this paper, we present observational results of molec¬ 
ular lines i 2C0(2-1), 13C0(2-1), and Ci®0(2-1). The to¬ 
tal gas mass (or Nhj ) is traced by extinction derived from 
NIR observations, independent from the molecular lines. 
The area covered in our observation is mainly den se gas 


with Ay > 3 mag. With the NICEST technique (|Lom- 

][», 

ag, wni 


bardi 


we are able to trace extinction up to Ay 


33 mag, which is mu ch deeper than achie ved in a simi¬ 
lar study of Orion by [Ripple et al.| (|2013|) w ho used the 
NICER technique (Lombardi & Alves||2U(Jl I to trace ex¬ 
tinction depths up to Ay ~ 15 mag. We also observed the 
optically thinner C^®0(2- 1) line, which provides insights 


into CO chemistry (e.g., Lada et al. 1994 Shimajiri et 
al.|2014 1 that cannot be readily derived usirig only CO 


and ^®CO observations. 

In the following, we describe our observations and data 
reduction in Section!^ We present the results and anal¬ 
ysis in Section Discussions and conclusions will be 
presented in Section]^ and Section]^ respectively. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. Near-Infrared Wide-Field Mapping Survey 
2.1.1. Telescope and Observations 

Near-IR observations for this study were obtained at 
the 3.5 m telescope of the Centro Astronomico Hispano 
Aleman observatory (hereafter CAHA) at Calar Alto in 
Almerfa, Spain. Our observations were made with the 
OMEGA 2000 camera, which has a field of view of 15'. 
In this paper we report observations for eight fields that 
cover the L1482 region. These fields are part of a more 
comprehensive survey of the GMG with a total of 24 
fields. Fields CNOl to GN08 (Table [B were observed 
in the three main near-IR broadbands, j, H and ATg, be¬ 
tween October 2009 and January 2010, with acceptable 
weather and seeing, except for one held, GN03, which 
was repeated on December, 2012. All observations con¬ 
sisted of 20 exposures of 60 seconds, coadded as 1200 
second exposure mosaics. The helds observed with the 
GAHA 3.5 m telescope were selected to cover the main 
hlament or “spine” of L1482, including the region of the 
LkHa 101 cluster. The pixel scale of OMEGA 2000 at 
the 3.5 m telescope is 0.45 arcsec pix“^, with excellent 
uniformity and negligible geometric distortion across the 
FOV. 

A list of all helds observed and considered for hnal anal¬ 
ysis can be consulted in tablej^ The table lists the held 
identihcation, the center of field positions, observation 
date, hlter, an estimate of the seeing based on the aver¬ 
age FWHM of the stars detected in each held, and the 
peak values for the brightness distributions, which are a 
good measurement of the sensitivity limits achieved. Be¬ 
low, we describe the observations and the data reduction 
process, including the construction of the photometric 
catalogs used to make our dust extinction map. 


2.1.2. Near-Infrared Data Reduction and Calibration 

The near-IR image data from OMEGA 2000 was re¬ 
duced with pipelines that made use of standard IRAF 
procedures complemented with stand alone routines. 
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Table 1 

Near-Infrared Observations of L1482 Fields 


Field ID 

Date Obs. 

Center Coords. 

J2000 

Filter 

Seeing 

[(")] 

LF Peak'^ 
[mag] 

(JALAH AL'i'O 3.5 m UMEGA 2U0U UBJSEHVAi'iOJNIS 

CN-Ol 

2009-12-24 

04:30:58.81 

-1-34:51:57.5 

J 

1.49 

19.00 

CN-01 

2009-12-24 

04:30:58.81 

-1-34:51:57.5 

H 

1.52 

19.75 

CN-Ol 

2009-12-24 

04:30:58.81 

-1-34:51:57.5 

Ks 

1.39 

19.25 

CN-02 

2009-12-08 

04:30:38.55 

-1-35:04:20.4 

J 

1.36 

19.50 

CN-02 

2009-10-09 

04:30:38.55 

-1-35:04:20.4 

H 

1.59 

20.75 

CN-02 

2009-10-09 

04:30:38.55 

-1-35:04:20.4 

Ks 

1.47 

19.75 

CN-03 

2012-12-23 

04:30:07.84 

-1-35:16:06.6 

J 

1.11 

22.25 

CN-03 

2012-12-23 

04:30:07.84 

-1-35:16:06.6 

H 

1.06 

20.25 

CN-03 

2012-12-23 

04:30:07.84 

-1-35:16:06.6 

Ks 

1.12 

19.25 

CN-04 

2009-10-07 

04:30:44.04 

-1-35:27:49.5 

J 

1.05 

21.25 

CN-04 

2009-10-07 

04:30:44.04 

-1-35:27:49.5 

H 

0.96 

19.50 

CN-04 

2009-10-09 

04:30:44.04 

-1-35:27:49.5 

Ks 

1.05 

21.50 

CN-05 

2009-10-07 

04:30:44.04 

-1-35:41:24.1 

J 

0.91 

21.25 

CN-05 

2009-10-07 

04:30:44.04 

-1-35:41:24.1 

H 

0.93 

22.00 

CN-05 

2009-10-06 

04:30:44.04 

-1-35:41:24.1 

Ks 

1.17 

19.75 

CN-06 

2009-10-06 

04:30:44.73 

-1-35:54:54.8 

J 

1.29 

22.00 

CN-06 

2009-10-06 

04:30:44.73 

-1-35:54:54.8 

H 

1.26 

21.00 

CN-06 

2009-10-06 

04:30:44.73 

-1-35:54:54.8 

Ks 

1.30 

20.00 

CN-07 

2010-01-02 

04:29:08.58 

-1-36:29:18.7 

J 

1.07 

20.75 

CN-07 

2010-01-02 

04:29:08.58 

-1-36:29:18.7 

H 

1.13 

19.00 

CN-07 

2010-01-02 

04:29:08.58 

-1-36:29:18.7 

Ks 

1.18 

18.75 

CN-08 

2010-01-03 

04:28:16.37 

-1-36:29:58.4 

J 

1.37 

20.25 

CN-08 

2010-01-03 

04:28:16.37 

-1-36:29:58.4 

H 

1.35 

20.00 

CN-08 

2010-01-03 

04:28:16.37 

-1-36:29:58.4 

Ks 

1.02 

19.25 


^ Expresses turnover point of observed magnitude distribution 


These pipelines are based on the FLAMINGOS near- 


Roman-Ziihiga 

2006 

Levine 

2006 


is essentially ident ical to the one described in [Roman- 1 


Zuniga et al. (20101 and we refer the reader to that paper 
tor details. 

Both photometry and astrometry of OMEGA 2000 
data products were calibrated relative to Two-Micron 
All Sky Survey (2MASS) lists obtained from the All-Sky 
Point Source Catalog (PSC). The final photometry cata- 
logs were mer ged with a combination of TOPCAT-STIL 
( Taylor|[2005 1 and IDL routines. In the case of adjacent 
helds, overlapping areas were treated with a routine de¬ 
signed to handle duplicate detections in such a way that 
we list preferentially a higher quality observation (e.g. 
a smaller photometric uncertainty) over a lower quality 
one. 

The 2MASS data were also used to replace observa¬ 
tions for saturated stars in all our frames. We used, in 
all cases, lists retrieved from the 2MASS All-Sky Point 
Source Catalog (PSC) at the Infrared Processing and 
Analysis Center (IPAC). The final photometry catalogs, 
containing either J, H, and Ks, or H and Ks photome- 
try were mer ged with a combination of TOPCAT-STIL 
( TaylorpOOS ) and self-made IDL routines. 


2.2. The CO Molecular-Line Survey 
2.2.1. Telescope and Observations 

A series of observations were carried out with the Hein¬ 
rich Hertz Sub-millimeter Telescope (SMT) on Mount 
Graham, Arizona, from November 2011 to April 2012. 
The SMT is located at an elevation of 3200 meters. The 
receiver has prototype ALMA Band 6 sideband separat¬ 
ing mixers with two orthogonal polarizations. Typical 
single sideband system temperatures during the observ¬ 
ing were around 200 K. The molecular lines ^^CO(2-l) 


(230.538 GHz), i3CO(2-l) (220.399 GHz), and C^^O{2- 
1) (219.560 GHz) were observed over a selection of 
17 10' X 10' square regions (hereafter “tiles”) along 
the south-eastern dense ridge of the CMC, roughly at 
RA=04 h30"^00L DEC=35°50'00" (see Figurealso fig¬ 
ure 1 of Lada, Lombardi, fc Alves|2009 ). Thes^iles have 
Ay ^ 3 mag, and cover the higher extinction regions of 
the cloud as well as the LkHa 101 cluster (roughly cen¬ 
tered in tile 12). 

In each tile, we performed the observations using the 
“on-the-fly (OTF)” mode. The scanning rate was 10" 
per second along the direction of Right Ascension (RA), 
with each row being separated by 10" in Declination 
(DEC). Each sampling point has a 0.4-second integra¬ 
tion, corresponding to 4" spatial sampling in RA with 
a beam size around 35", therefore the source is fully 
sampled. The total integration time for each tile is 
roughly 1 hour. We utilized dual side band mode and 
had two distinct sets of observations: Obsl: ^^CO(2- 
1) in lower side-band (LSB) and ^^CO(2-l) in upper 
side-band (USB); Obs2: Ci®0(2-1) in LSB and 1^00(2- 
1) in USB. In each set, two orthogonal linear polariza¬ 
tions were observed simultaneously. The total band¬ 
width we cover with the spectrometer is 128 channels 
of 0.25 MHz each or 32 MHz total, for each line (^^CO, 
^®CO, and C^®0). The channel resolution is 0.25 MHz, 
set by the filter widths, and this corresponds to about 
0.34 km s“^ velocity resolution, depending on the line 
rest frequencies via the doppler formula. We checked 
the pointing on the compact nearby CO source CRL618 
(RA=04*'43“42L9, DEC=36°08'14.7") at the beginning 
of every tile. Five-second integrations were performed 
on each point of a 5-point cross pattern centered on that 
source. The pointing error was typically < 5". To cal¬ 
ibrate the intensity, we observed the strong CO source 
W30H (RA=02*'27“59L2, DEC=61°55'43.2") between 
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Figure 1. Contour map of visual extinction derived from our 
near-infrared survey with the NICEST method. Contours at inter¬ 
vals of Ay = ( 5, 7, 10, 12, 15, 20, 25 and 30 mag). See descriptions 
in Section |2.1| The square boxes outline the 17 tiles representing 
the areas surveyed in CO emission. LkHa 101 is marked by the 
star in tile 12. 

every two tiles. 

2.2.2. Millimeter-wave Data Reduction 

All data were first processed through the CLASS soft¬ 
ware package. We fitted linear baselines and subtracted 
them from the data. The data image cubes wer e ported 
into MIRIAD software format (|Sault et al.|1995 ) for fur¬ 
ther analysis and display. OTF data were interpolated 
onto regular grids with grid spacing of 10" before being 
stored into the MIRIAD data set. Next the data were 
calibrated using W30H. A scale factor is derived through 
the observation toward W30H between every two tiles. 
The factor is used to scale the antenna temperatu re T1 
to the main-bea m brightness temperature T„ib (see|Bieg~ 


ing et al. 2010 for details). For each molecular tran¬ 
sition, the two polarizations were averaged to increase 
the SNR. Finally all 17 tiles were combined together to 
form a single map (weighted by their RMS noise at the 


small overlapped area). The original data cube, with ^ 
0.34 km s“^ velocity channels (set by the 0.25 MHz filter 
bandwidth) were re-sampled onto a 0.15 km s“^ spacing, 
by 3rd-order polynomial interpolation. This re-sampling 
puts all of the CO isotopologues on exactly the same 
velocity grid, so we can make ratio maps, for example. 
We convolved all maps with Gaussian kernels (FWHM 
= 18.1" for map, FWHM = 13.6" for i^CO and 
C^®0 maps) to have a final map resolution (FWHM) of 
38" with a grid spacing of 19" in order to match the res¬ 
olution of and facilitate direct comparison with our dust 
extinction maps. 

We calculated the RMS noise for all emission-free chan¬ 
nels in each map to be about 0.1 K for all three lines. 
We compared the integrated intensity maps of the two 
sets of ^^CO maps by subtracting one from the other. 
We found a small spatial difference (typically < 10" be¬ 
tween each pair of corresponding tiles), and we shifted 
and re-gridded Obsl ^^CO and ^^CO maps to Obs2. We 
averaged the two sets of ^^CO maps weighted by their 
RMS noise to lower the noise in the averaged ^^CO map 
to 0.07 K (RMS). 

3. RESULTS AND ANALYSIS 
3.1. The NICEST Dust Extinction Map 

Figure shows the deep dust extinction map of the 
southeast portion of the CMC derived from our observa¬ 
tions. For this map Calar Alto data were supplemented 
with data from 2MASS in some outer portions of the 
surveyed area where Calar Alto data were not available. 
Also shown are the boundaries of the 17 tiles of the 
multi-line CO mapping survey. The extinction map was 
constructed with the optimized version of the Near In¬ 
frared Color Excess Revised (NICER) technique, known 
as NICEST (Lombardi & Alves 2001 Lombardi] 2009). 
The NICER/JMCES'L technique allows us to measure 
dust extinction from the infrared excess in the colors of 
background stars caused by dust absorption. The excess 
is derived by assuming as intrinsic colors, the colors of 
stars detected in a nearby off-cloud control field with neg¬ 
ligible extinction (see Table 1). Eor the measurements of 
extinction, we avoid the use of sources with intrinsic color 
excess, such as candidate young stellar objects and dusty 
galaxies. NICEST is essentially identical to NICER, ex¬ 
cept that the Ay estimator is modified to account for 
small-scale inhomogeneities, due to a bias introduced by 
the decrease in the number of observable stars as ex¬ 
tinction increases toward the denser areas of a molecular 
cloud. The bias is removed by modifying the estimator 
with an additional term that accounts for the expected 
decrement. 

Extinction measurements for individual sources are 
smoothed with Gaussian filter with a FWHM of 38", and 
the maps are constructed with 19" spatial sampling. At 
each (line-of-sight) position in the map, each star falling 
within the beam is given a weight calculated from a Gaus¬ 
sian function, and the inverse of the variance squared, 
which is in turn also used for the bias-correction factor. 
A preliminary value of Ay at the map position is calcu¬ 
lated as the weighted median of all possible values within 
the beam. Then, we estimate the mean absolute devia¬ 
tion (m.a.d.) of these values to remove large deviates 
using m.a.d.-clipping. The remaining values are used by 
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the NICEST estimator to calculate the final Ay at the 
corresponding position. 

In the case of pixel positions where only 2MASS 
sources are available, the number of sources per beam 
is significantly smaller, and there are cases where there 
are not enough stars to estimate Ay; in those cases, our 
code increases the size of the beam by a factor of 2 over 
the nominal value. The nominal 38" FWHM of our maps 
represents an increase of resol ution by a factor of 2.1 com¬ 
pared to the earlier study of Lada, Lombardi, &: Alves] 
(2009). This resolution is comparable to the beam size 
of the CO maps. The resulting noise of extinction map 
is 0.05 mag for Ay < 10 mag and 0.1 mag for Ay > 10 
mag. 


3.2. CO Spectra and Integrated Intensity Maps 

In order to provide an overview of the observations and 
the general quality of the data we first show in Figure 
the averaged CO spectra derived for each of the 10' x 10" 
tiles of our mapping grid. On this coarse scale (1.3 x 
1.3 pc) the ^^CO profiles exhibit evidence for both self¬ 
absorption and multi-component structure. The pres¬ 
ence of self-absorption indicates that the ^^CO emission 
is quite optically thick. The profiles of the rarer iso¬ 
topes appear to be less complex, mostly characterized 
by simple Gaussian-like shapes. Maps of the emission 
from 12C0(2-1), i3cO(2-l), Ci®0(2-1) integrated over 
velocity ranges, —12 to -1-10 km s“^, —10 to -1-6 km s“^, 
—4 to -1-3 km s“^, respectively, are shown in Figureas 
grey scale maps with contours overlaid for clarity. The 
maps display some common morphological features. The 
strongest emission occurs in the mid-southern area (tiles 
11-14, 6 « 35.25°) and is in close proximity to the LkHa 
101 cluster (which lies within tiles 10-12, S « 35.27°). 
On close inspection of the maps, a sharp drop off in CO 
emission is observed along a warm Hlamentary structure 
southwest of the cluster. This sharp drop in CO emis¬ 
sion together with the enhanced CO emission within the 
warm filament likely indicates a physical interaction be¬ 
tween the cluster stars and the filamentary ridge struc¬ 
ture. A similar structure with an accompanying drop off 
is also present in the extinction map of Figure 1. Here 
we can see that the warm filamentary structure is a small 
section of the more extended high extinction ridge that 
forms the backbone of the cloud in this region and which 
is the subject of this survey. In these more extended 
regions away from the cluster, the CO emission is com¬ 
paratively quiescent. 

The RMS noise for the integrated intensity maps aw 
(in K km s“^) was first estimated using 

aw=acX ( 1 ) 

V 'channel 

where Uc (in K) is the RMS noise per channel derived 
for emission-free velocities, AV is the velocity width over 
which the integrated intensity is calculated. Nchannei here 
is the number of real independent channels or filters (of 
bandwidth 0.25 MHz, ~ 0.34 km s“^) within AH. The 
estimated aw for ^^CO(2-l), ^^CO(2-l), C^®0(2-1) in¬ 
tegrated intensity maps are 0.16, 0.22, 0.16 K km s“^, 
respectively. We also calculated the RMS for an inte¬ 
grated intensity map computed over the same number of 
emission free channels in the data cube. The aw derived 


this way are 0.11, 0.24, 0.24 K km s respectively. To 
be conservative we adopt the larger value for aw, i-e., 
0.16, 0.24, 0.24 K km s-\ for 12C0(2-1), i3cO(2-l), 
C^®0(2-1), respectively (hereafter aw,12, o'w,i3, crvF.is)- 
The ^^CO and ^^CO lines are very strong and in our 
subsequent analysis we only consider integrated intensi¬ 
ties detected above baw- With C^®0, we only consider 
detections above 2aw- We also note some striping ap¬ 
parent in tiles 5, 16 and 17 of the C^®0 map. These 
are likely artifacts due to system temperature variations 
in some of the OTF scans in those regions of the cloud. 
Prior to the subsequent analysis of these data reported 
below this map was “cleaned” so that pixels containing 
the stripes were removed from the C^®0 databas^ 

3.3. CO excitation and column densities of CO and 

C^^O 

In this paper we adopt the standard LTE analysis to 
determine excitation temperatures and column densities. 
The main beam brightness temperature, Tmb) is given by 


Tmb = [MTe.) - MTbgm " e”"") 


( 2 ) 


where Jv{T) = is the Rayleigh-Jeans equiva¬ 

lent temperature of a blackbody of physical temperature 
T, Tex the excitation temperature of the line under con¬ 
sideration, Tbg the temperature of the background (in 
this case the CMB), and the optical depth of the tran¬ 
sition. The excitation temperature Tex of a given line can 
be derived from the peak temperature of the line, pro¬ 
vided the line is optically thick. Typically ^^CO emission 
is optically thick in GMCs. In the CMC a large opacity 
for the ^^CO(2-l) line is confirmed by the ^^CO/^^CO 
integrated intensity ratio. The probability density distri¬ 
bution for this ratio for all pixels is found to peak near 
a value of 3.5, considerably s maller than the optically 
thin limit of « 40-70 (see, e.g.,|Langer & Penzias||1993|). 
Less than 10% of the pixels are characterized by ratios 
exceeding 10. This is consistent with the fact that the 
surveyed area was dominated by relatively high extinc¬ 
tion (dust column density). For t 1, Tex for the ^^CO 
line is calculated from the following equation (by solving 
Eq§: 


T — 

-J- per. — 


hv/k 


ln{l + 


hv jk 

{Tbg) 


(3) 


where T„ib,p is the peak brightness temperature in the 
^^CO line. At this frequency, J,,{Tbg) = 0.19 K and 
is generally much less than the peak line temperature 
(Tmb,p)- The typical erroij^in Tex propagated from the 
3 sigma uncertainty in Tmb,p is 0.22 K. The maxi¬ 
mum percentage error in Tgx is below 5%. Note that in 
some regions, e.g., tiles 12 and 14, the ^^CO line shows 
self-absorption, this would result in underestimation of 
Tmb,p and thus Tex- For Tnib,p between 6 K and 50 K, a 
change of this quantity by 1 K produces a corresponding 
change of Tex by 1 K. 


^ In all cases, the stripes show up in only one polarization. We 
masked them before combining the two polarization maps. The 
striping area would therefore have higher noise (by a factor of 1.4). 
^ This value encloses 90% of pixel errors. 
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Figure 2. Averaged spectra for all three molecular lines in each 10' X 10' tile. Roughly 1000 spectra (from 1000 pixels) in each tile were 
averaged together. C^®0{2-1) spectra were enlarged by a factor of 3 to be shown together with the other two line spectra. 


Figure shows a map of the spatial distribution of 
the excitation temperatures in the surveyed region. 
Over the vast majority of the cloud, Tex is found to be in 
a narrow range around 10 K as expected for gas heated 


by cosmic rays and cooled by molecular li nes (Goldsmith 


& Langer 1978 


Bergin & Tafalla 20071. Higher tem¬ 
peratures are found in the southern part of the region. 
The highest excitation temperatures approach 40 K and 
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Figure 3. Integrated intensity maps of (a) ^^CO(2-l), (b) '^®CO(2-l), (c) C'^®0(2-1) of the Southeast region of the cloud shown in grey 
scale with contours overlaid for clarity. The integrated velocity ranges are -12 to 10 km s“^, -10 to 6 km s“^, -4 to 3 km s“^, respectively. 
The maximum integrated intensities are 130.0, 51.5, 7.1 K km s“^ for (a) ^®CO(2-l), (b) ^®CO(2-l), (c) C^®0(2-1), respectively. The 
contour levels are 10%, 20%, 40%, 60%, 80% of the maximum of the respective maps. All the maps have been placed on the same spatial 
grid with a resolution of 38” sampled at 19”. See text for more details. 


are found within the bright filamentary ridge southwest 
of the LkHa 101 cluster (see figure) suggesting that the 
cluster is providing additional heating of the gas in that 


area. 


The optical depths of ^^CO(2-l) (ris) and C^®0(2-1) 
(ris) were calculated assuming LTE from: 


= —ln\l — 


T 

-J- 7 


mb 


hi'/k 


l)-i - - 1) 


-1 


(4) 


Here Tex for both transitions was taken to be that de¬ 
rived for ^^CO(2-l). The probability density distribution 
of opacities for ^^CO peaks around 0.5 with the majority 
(87%) of pixels appearing to be relatively optically thin 
(i.e., T < 1). For C^®0 the distribution peaks around 
0.2, significantly less than that for ^^CO (as expected), 
and with r < 1 for essentially all pi xels in the map. 

We followed Pineda et al. (2010) in calculating both 
the ^^CO and C^®0 column density using Eqj^for J=2-l 
transitior0 The rotational constant Bq and the Einstein 
coefficient A. 2^1 are taken from an online databas^ The 
partition function Q is estimated by EqU which is ac¬ 
curate to 10 percent for T > 5 K (note that in our data. 


Tex > 5 K everywhere, as shown in Figure]^. 


N = 




'^bg — 2 


Q 


L ^3 A ^ jy —h Bq J (^J1) 

hc^A 2 ^i (2 J + l)e—- 

J T{v)dv f ^ ^ 

/[l_e-HD]dW ) 


-/iSn J(J+1) 

Q = Ej(2J+l)e 


hBn 


(5) 


( 6 ) 


We note that sub-thermal excitation for the optically 
thin ^^CO(2-l) and C^®0(2-1) lines might be expected to 
be a concern in low extinction regions of the cloud where 
number densities may be lower than the critical density 
for thermal (LTE) excitation (« 3 x 10^ cm“^). This 
could lead to underestimates of the CO column densi¬ 
ties in the outer regions. However, since the studied area 
here is dominated by relatively high extinction (Ay ^ 
3-4 magnitudes) we expect that the effects of any sub- 
thermal excitation are minimized. We can roughly esti¬ 
mate the magnitude of this effe ct using the RADEX, non - 
LTE, radiative transfer code (van der Tak et al. 2007). 


We find that with n(H 2 ) = 5 x 10^ cm Tgas = 10 K 
and W(i3co) « 2.5 K km s-\ Nradex « 2.0 Nlte- 


^ We noticed t he typo in equation (18) iniPineda et al.| ||2010||. 4, DISCUSSION 

as pointed out by [Ripple et al.|ll2013ll 

^ http://home. strw. leidenuniv.ni/~inoida’ta/co.htmi 4 . 1 . CO ColuTTiTi Densities and \^isu(il Extinction 
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Figure 4. Tex map of the Southeast region of the cloud derived 
from CO(2-l) map using Equation]^ The boxes show the observed 
tiles 5-17. The blue open circle shows the position of the LkHo 101 
cluster (from SIMBAD). 


4.1.1. CO and Abundances 

Since our measurements of CO column densities and 
dust extinction are set to the same angular resolution 
and spatial grid, we are in a position to directly compare 
the two data sets. We note here that the angular scale of 
a pixel is 19", corresponding to a spati al scale of ^ 0.04 
pc 8000 AU ) at a distance of 450 pc (Lada, Lombardi, 
fc Alves||200^ . 

Assuming a constant gas-to-dust mass ratio. 


(N ( H)/Av =1.88x10^^ cm~^ mag~^) (iBohlin e t 
ah] |1978| iRachford et al.| j2009| |Pineda et al.||2UU8|), 
the CO abundance can be expressed as: [CD] = 
N(CO)/N(H) = N(CO)/[Av(1.88 x lO^i cm-^ mag-^)]. 
For optically thin emission lines, such as and 

C^®0, we expect the gas column density to increase 
linearly with Ay, provided the clou d is characterized 
by a constant CO abundance_(e.g., |Lada et al.||1994 
Alves et al. 1999). In Figure we plot the molecular 
column density vs. extinction at each pixel in the 
mapped region for ^^CO (a) and C^®0 (b), respectively. 


Contrary to our simple expectations, each plot exhibits 
considerable scatter with only weak correlations between 
gas and dust column densities. If the gas in the CMC 
were characterized by a single, constant abundance, 
the relation between its column density and extinction 
would consist of a single straight line in each plot. At 
fixed Ay there would be only small scatter in N(CO), 
which is clearly not the case. For instance, at Ay ^ 
20 mag the ratio between the highest and lowest values 
of N(^^CO) and N(C^®0) is > 10 for each isotope. 
Clearly, the two plots suggest strong spatial variations 
in abundances of these two species across the cloud. 

On closer inspection of Figure one can almost make 
out two signihcant branches (correlations) of pixels in 
both plots. First, a large group of pixels in each plot 
(below the blue dashed line, i.e., N(^^CO) < 1.0x10^® 
cm“^, N(C^®0) < 1.0x10^® cm“^) display an almost 
constant, low column density across a wide range of Ay. 
This flat relation between the two quantities indicates 
a trend of decreasing molecular abundance ([^^CO] and 
[C^®0]) from low to high Ay pixels. Second, another 
group (above the blue dashed line) of pixels is character¬ 
ized by gas column densities that appear to rise more or 
less linearly with extinction to very high column density 
at Ay ~ 20 mag (I.OxlQi® cm-^ < N(i3CO) < e.OxlQi® 
cm“^, l.OxIO^® cm“^ < N(C^®0) < 6.0x10^^ cm“^), 
indicating a more or less uniform molecular abundance 
with depth into the cloud for this group of pixels. These 
branches likely mark two limiting extremes of spatially 
dependent variation in the chemical conditions within 
the CMC. 


Spatial variations in ^^CO abundances have been re¬ 
ported previously in a few other molecular clouds, al¬ 
though in these clouds the variations appear less pro- 
nounc ed than found here. For the Perseus molecular 
cloud, Pineda et al. (2008) found generally tighter cor- 
relations of ^'^CO column density with extinction than 
found here, although they explored the correlation at 
much lower extinctions (i.e.. Ay < 10 mag) than in this 
study. Nonetheless, over this lower extinction range they 
did find measurable variations (at the ~ 40% level) in 
the abundances with posi tion across the clou d. For the 
Orion A molecular cloud. Ripple et al. (2013) measured 
abundances at both low and high extinction and found 
more signihcant positional variations in the ^^C O abun¬ 


dance acros s the cloud than observed in Perseus (Pineda 
et al. 2008). At Ay ~ 20 mag, N(^^CO) in Orion A 
was found to vary by about a facto r of 4 between vari- 
ous locations in the cloud. Further, Ripple et al. (2013) 
were able to associate these positional variations in abun¬ 
dance with positional variations in the physical condi¬ 
tions within the cloud. At the lowest extinctions (Ay < 
3 mag), where there is not enough ^^CO to self-shield 
against UV dissociation, they found extremely low 
abundances with signihcant (factor of 8) spatial varia¬ 
tions. At intermediate extinctions (3 < Ay < 10 mag), 
where self-shielding is considerably more effective, higher 
abundances were found with modest (factor of 2) spatial 
variations. At the highest extinctions they observed the 
highest abundances but also the largest spatial variations 
in abundances across the cloud. This latter regime in¬ 
cluded cold portions of the cloud where CO depletion 
depressed the abundances as well as hot regions heated 
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Figure 5. A pixel-by-pixel comparison between (a) ^^CO, (b) column density and visual extinction Ay for the entire surveyed 

area. Each point represents a 19^^ pixel which corresponds to ~ 0.04 pc at a distance of 450 pc. A 6 (Jiv ,13 cut is applied to W(^^CO) in 
N(^^CO) calculation, while 2a\Y,i3 cut is applied to W(C^^O) in N(C^®0) calculation The blue dashed lines tentatively divide the points 
into two groups: an approximately constant abundance group above the line and a constant column density group below. The size of the 
erro r bar enclose s 90% of the calculated 3 sigma errors in column density (see text). The error of Ay is too small to be shown (Section 
See Section |4.1.1| 


by young stars where CO desorption produced enhanced 
CO abundances. 


4.1.2. Mapping Spatial Variations in CO Abundances 

The molecular abundance is regulated by the physi¬ 
cal and chemical conditions in the CMC. For instance, 
the gas-phase molecules can condense onto dust grains 
to form ice mantles, and the molecules on grain surfaces 
can be returned to the gas via thermal evaporat ion (e.g., 
van Dishoeck et al.|ll993 Caselli et al. 1999). These 
processes are tightly linked to the local temperature (es¬ 
pecially the dust temperature), and are responsible for 
setting the gas-phase molecular abundance. At low ex¬ 
tinctions where column densities of dust and gas are low, 
abundances are very sensitive to FUV radiation which 
can drive volatile chemistry via processes of fractionation 


and selective photo dissociation ( [Lada et al.||1994jjRollig 
fc Ossenkopf||2013). The large scatter in h'igure could 


plausibly be caused by strong variation of such physi¬ 
cal/chemical conditions in different sub regions of the 
surveyed area of the cloud. If this is the case, we would 
expect to observe a tighter correlation between column 
density and Ay if we restricted the plots to cover smaller 
spatial areas within the surveyed region, since the physi¬ 
cal conditions in such smalle r regions wou l d be e xpected 
to be more uniform. Indeed, Ripple et al. (2013) demon¬ 
strated the connection between such physical conditions 
and spatial CO abundance variations in the Orion cloud. 

In order to examine the distributions of N(^^CO) and 
N(C^®0) versus Ay on smaller spatial scales, we con¬ 
structed these relations for each individually observed 
10' X 10' mapping tile (c.f. Figures 1 and 2). In Figures]^ 
andmwe plot the results for ^^CO and C^®0, respectively. 
In the individual plots contained within these figures we 
indeed see a significant reduction in scatter compared to 


Table 2 

Survey Tile Properties 


Tile 

Tex 

(K) 

slope®' 

(10^® cm“^ mag“^) 

r-value®’ 

01 

9.21 

0.014 

0.42 

02 

9.31 

0.027 

0.65 

03 

8.30 

0.031 

0.64 

04 

8.52 

0.025 

0.50 

05 

9.14 

0.024 

0.57 

06 

9.23 

0.018 

0.43 

07 

9.99 

0.046 

0.64 

08 

9.19 

0.037 

0.55 

09 

12.3 

0.034 

0.37 

10 

12.3 

0.066 

0.67 

11 

13.1 

0.148 

0.80 

12 

19.6 

0.193 

0.76 

13 

16.5 

0.204 

0.91 

14 

16.9 

0.102 

0.62 

15 

9.07 

0.059 

0.75 

16 

8.65 

0.054 

0.64 

17 

9.66 

0.022 

0.64 


^ Linear regression results from Figure In tiles 
11-14, the fitting is performed over the entire Ay 
range. In the rest of tiles, the fitting is up to Ay 
= 10 mag. The r-value (Pearson coefficient) is an 
estimation of correlation coefficient. 

that characterizing Figure More importantly, there 
is very strong spatial variation in the relation between 
extinction and CO column density across the surveyed 
region of the CMC. This is evident in Figure 6 where we 
plot linear fits to the data. The results of the fits are 
shown in Table 2 where the slopes vary by nearly an or¬ 
der of magnitude across the cloud. This is also readily 
apparent in Figure where the axes on the plots have 
the same scaling and especially clear on close inspection 
of Figure]^ where the scale of the vertical axis ranges by 
as much as an order of magnitude between the various 
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Figure 6. A pixel-by-pixel comparison between column density and Ay within each tile. The scale of the column density axis 

varies between tile rows and a fiducial horizontal line representing a constant N(^®CO) of 0.5 X 10^® cm“^ is drawn for inter comparison. 
The individual points are color-coded by Tex. A 5 cr\Y,i3 cut is applied to W(^®CO) in N(^®CO) calculation. The red solid lines show linear 
regression in Ay > 3 mag for tiles 11-14, and in 3 mag < Ay < 10 mag for the rest. As a reference, the fitting result for tile 13 is shown 
as green lines in every tile. The magnitude of the error is the same as in Figure]^ a). 


tiles. These spatial variations in N(^^CO) and N(C^®0) 
give rise to the large dispersion in these quantities in the 
merged relation constructed for the entire observed re¬ 
gion as in Figure Despite the variation in the range of 
column densities at different locations across the cloud, 
the range in extinction spanned by the observations in all 
the plots is very nearly the same. The large variation of 
column density of and C^®0 with dust extinction 


confirms that the abundances of the two isotopologues 
exhibit large spatial variations. Consider Figure The 
tiles can be placed into two groups, one containing tiles 
1-10 and 15-17, and the other consisting of tiles 11-14. 
In the former group, N(^^CO) rarely exceeds a 5-6 xlO^^ 
cm“^, while in the latter group N(^^CO) can reach val¬ 
ues as high as 3-5 x 10^® cm“^ over the same extinction 
range. The N(C^®0) plots in Figure exhibit similar 
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Figure 7. A pixel-by-pixel comparison between column density and Ay within each tile. A 2(Tiy ig cut is applied to the integrated 

intensity in N(C^®0) calculation. The magnitude of the error is the same as in Figure[^b). The axes have the same horizontal and vertical 
scales for all tiles, otherwise the same as Figure]^ 


behavior. 

4.1.3. The Relation between Gas Temperature and CO 
Abundances 

From comparison of Figures [^and[^ with Figure H we 
see that the tiles (11-14) with^e largest values ofCO 
column density spatially coincide with the region of ele¬ 
vated ^^CO excitation temperature. Since the ^^CO(2-l) 
line is generally so optically thick, and since our sampled 
region preferentially contains high column density (Ay 


> 3-5 mag) gas, sub-thermal excitation of ^^CO is un¬ 
likely. Therefore the excitation temperature map traces 
the gas kinetic temperature and the distribution of the 
molecular abundances [^^CO] and [C^®0] are tightly re¬ 
lated to local physical conditions, in particular the gas 
kinetic temperature. As a further test of the connection 
between column density and temperature we have color- 
coded the points in Figures^ and [^according to the gas 
temperature in each individual pixel. Inspection of the 
figures shows that the pixels with the highest gas temper- 
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atures (Tgx ^ 15 K) are also the pixels with the largest 
CO column densities. Pixels with low gas temperatures 
(Tex ^ 15 K) rarely reach high CO column densities, in¬ 
dependent of extinction. These correlations imply that 
the abundances of ^^CO and C ^®0 are physically related 
to the gas temperature distribution. 

To examine more directly the relation between the 
molecular abundances and the gas temperature, we plot 
the log([^^CO]) - log(Tex) relation and log([C^®0]) - 
log(Tex) relation in Figure Both molecular abun¬ 
dances show dependence on Tex, although [^^CO] is ap¬ 
parently sensitive to Tex over a larger range of tempera¬ 
ture than [C^®0]. This confirms that [^^CO] and [C^®0] 
are both regulated by gas temperature, but over different 
ranges of temperature in the extinction range (3 mag ^ 
^ 25 mag) we observed. To quantify the relations, 
we empirically fit a power-law functions to the unbinned 
data for both correlations. Tex is derived from the ^^CO 
line, which has very good SNR. Since the noise of ex¬ 
tinction map is relatively small, i.e., 0.05 mag for Av<10 
mag and 0.1 mag for Av>10 mag, uncertainties in the 
abundances are dominated by the noise in the molecular 
line data. 

The fits were performed over different ranges of Tex 
for the two lines: Tex > 10 K for ^^CO and Tex > 17 K 
for C^®0. The low temperature cutoffs were empirically 
chosen, by inspection, to correspond to the ranges where 
the relations appeared to be linear. However, it did not 
escape our attention that the low temperature cutoff for 

/-»i-\T'T'z:iC'TA/~kn/^lo on Vvl in-f i/~vvi +-^TVTt-\£if n-f n t'Oi r\f- 


lAlves et al. 

1999 

Bergin et al. 

MIT 


Alves et al.||1998 2001 Bergin et al.||200r ) and the pres- 


ence of emission from high dipole molecules such as NH 3 
from such regions indicate that they likely correspond 
to sufficiently high gas volume densities (> 10 ^ cm“^) to 
promote CO condensation onto grain surfaces. 

In Figurej^we plot the relation between (a) and 

(b) C^®0 abundance and extinction in the CMC. Al¬ 
though there is a large scatter in both relations, it is 
clear that the abundances of the two species decline with 
extinction in a manner consistent with similar abundance 
reductions seen in other molecular clouds and suggestive 
of CO freeze out onto grain surfaces in the interior re¬ 
gions of the cloud. If this decline in abundances at high 
extinctions is due to freeze out, we also expect the reduc¬ 
tion of the gas-phase abundances to be occurring in the 
coldest regions of the cloud. In Figure P the points are 
color-coded by the gas temperature and we see both a 
preponderance of cold gas in the cloud, and a prominent 
decrease in abundance with extinction for the cold pixels, 
while the warm pixels show no decrease. This supports 
the hypothesis that the decline in abundances with ex¬ 
tinction in cold regions is largely due to depletion. Of 
course, depletion requires low dust temperature and our 
analysis has been based on measurements of gas temper¬ 
ature. Ga s temperature genera lly tracks dust tempera¬ 


ture (e.g., Forbrich et al. 2014) but direct observations 
of dust teinperature would make for a more compelling 


- 17 K ( 

van Dishoeck & Black 1988 

), and as discussed 

case. Fortunately, published maps of dust temperature 

later this 

may not be coincidental. 4' 

tie linear regression 

tor the CMC based on Herschel observations exist I 

Har- 


tion coefficient = 0.66; [C^®0] oc T^^®® with correlation 
coefficient = 0.34. We performed a test on the power- 
law model for [^^CO] - Tex and [C^®0] - Tex relations. 
The is calculated using: 


X = S 


(O - E)2 


(7) 


erally good agreement between the gas and dust tem¬ 
perature spatial distributions. In particular, in the cold 
regions away from LkHa 101 (corresponding to our tiles 
01-09 and 15-17, respectively), th e median Ta 14. 5 
K and the minimum Td ~ 10 K (Harvey et al. 2013), 
far below pure CO ice sublim atiori temperature 17 


K 


where O and E stand for observed and modeled [CO] 
(calculated from the fitted power-law), respectively, and 
Ua the local standard deviation derived in a 2 K bin en¬ 
closing the data. At a significance level of 0.05, the fitted 
power-law models for [^^CO] - Tex and [C^®0] - Tex re¬ 
lations are acceptable, and the x^ per degree of freedom 
are 1.02 and 1.06 for [^^CO] and [C^®0], respectively, 
meaning the fit is reasonable. The implications of these 
results will be discussed in the next two sections of the 
paper. 

4.1.4. Molecular Depletion and Desorption 

We argue here that the dependence of the ^^CO an^ 
C^®0 abundances on gas temperature shown in Eigure 
is a result of the combined effects of gas depletion am 
desorption in the CMC. Theoretical considerations sug¬ 
gest that sticking of CO onto dust grains can remove 
the molecule from gas-phase very efficiently within cold 
(T< 17 K) and dense (n(H^ > 10"* 


van Dishoeck e t al. 1993). Under such conditions, 
Bergin et al.| ( |l995 ) showed a relatively short CO deple¬ 
tion timescale, < 10® yr. So with Tj ~ 10 K and riHa = 
10"^ cm“®, CO would unavoidably deplete. This is con¬ 
sistent with the reduction of CO abundance we observe 
in the cold, high extinction regions in the CMC. 

At the higher gas temperatures (T > 17 K) we would 
expect that thermal evaporation (desorption) of the gas 
from the grain surfaces releases any trapped CO and in¬ 
creases its gas phase abundance relative to the cold re¬ 
gions of the cloud. Indeed, that is what is clearly shown 
in Figure]^ where the observed abundances increase with 
increasing gas temperature. The warmest CO gas in our 
surveyed region is found in four contiguous tiles (11-14) 
that are in close proximity to the massive star LkHa 101 
and Jts associated cluster (shown as the blue ci rcle in Fig¬ 
ure 
al. 


cm ® ) regions o f 


1993). 


molecular clouds (see, e.g., [van Dishoeck et al 
Indeed, compelling observational evidence tor reduction 
of [C^®0] has been reported for numerous cold (T < 10 
K) and high extinction (Ay ^ 10 magnitudes) regions of 


The CMC dust temperature map from Harvey et| 
(see their figure 4) shows a vast region of high 
Hust temperature (Td ^ 28 K) near LkHa 101, which 
coincides with the positions of tiles 11-14 in our obser¬ 
vation. It is in this warm area that the CO abundances 
appear to be at their highest. CO thermal evapor ation is 
caused by hi gh dust temperature. For instance, |Bergin| 
et al. (1995) reported that shortly after a “star turns 
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Figure 8. A pixel-by-pixel comparison between (a) ^^CO, (b) abundance and Tex for the entire surveyed area. Each point 

represents a 19^' pixel which corresponds to ~ 0.04 pc at a distance of 450 pc. Abundances were calculated from pixels where W(^^CO) > 
^^w,i3 and W(C^^O) > The yellow contours sh ow th e surface density of points. The blue lines show the linear fitting within Tex 

> 10 K for and Tex > 17 K for C^®0. See section [4.1.3| 
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Figure 9. The relation between CO abundance and extinction for (a) ^^CO, (b) over the entire surveyed area. Although there 

is large scatter, both plots show that the abundances of these two rare isotopologues clearly decrease with increasing extinction into the 
cloud. This is likely a result of CO depletion in the cold high extinction regions of the cloud. The individual points are color-coded by the 
value of the Tex derived from ^^CO observations at their position. Note that the hotter pixels have elevated abundances compared to the 
cold pixels and these abundances don’t appear to decrease with extinction. 

magnitudes, similar to those studied here for the CMC. 
They found that in regions with 5 < Ay ^ 10 magni- 
tudes, N(’^CO) was generally linearly correlated with 
Ay. However, in those partitions (2, 3 & 4) that also con¬ 
tained higher extinction material and were characterized 
by mean gas excitation temperatures that were below ~ 
22 K, the N(’^CO)-Av relation became flat above Ay 
« 10 magnitudes. On the other hand, in two partitions 
(5 & 6) with both high extinctions and high mean gas 
excitation temperatures (Tex ~ 26-33 K), N(’^CO) con¬ 
tinued to rise above Ay « 10 magnitudes. [Ripple et al. 
(20131 posited that this behavior of N(’^^CO) withextinc- 
tion was a result of depletion/desorption effects at Ay > 
10 magnitudes. Our res ults for the CMC ar e certainly 
consistent with those of Ripple et al. (2013) for Orion 
A. In this respect, the physical/chemical conditions in 
the Orion and the CMC appear to be quite similar. The 


on”, the dust temperature reaches Td = 25 K (riHa = 
10"^ cm“^, Tgas = 20 K), and all CO is in the gas phase 
within 100 yr. The abundance enhancement near LkHa 
101 is likely caused by hot dust h eated by the cluster 
Over most of the rest of the cloud Harvey et al. (2013) 
find Td 10 -14 K. In these regions the gas temperatures 
are also found to be near 10 K and a significant amount 
of CO should be locked in ice mantles on t he dust grains 


at lea st at the highest extinctions (also seefBergin et al. 
1995). These considerations suggest that the variation 
^CO] and [C’^®0] in the CMC is mostly regulated by 


of 


the chemistry of depletion and desorption which in turn 
is driven by the (dust) temperat ure distribution. 

Recently, [Ripple et al. (2013) have reported similar 
results for the Orion A rnolecular cloud. They divided 
that cloud into 9 spatial partitions (1-9) and considered 
regions within those partitions with extinctions. Ay > 5 
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spatial variations in abundances of the rarer isotopes of 
CO are largely determined by depletion/desorption ef¬ 
fects due to spatial variations in dust (and gas) tem¬ 
peratures resulting from localized star formation activity 
within the clouds. 


4.1.5. The ^^CO-to-C^®0 Abundance Ratio: Influence of 
UV photodissociation 

In the outer, lower extinction regions of molecular 
clouds, FUV rad iation is expected to play a ro l e in cloud 
chemistr y (e.g., van Dishoeck & Black 1988 Visser et 
al. 20091. In low extinction regions (i.e., Ay ^ 3 rnag) 
fractionation of can increase its abundance rela¬ 

tive to ^^CO. At somewhat higher extinctions (Ay ^ 5 
mag) selective photodissociation of C^®0 is likely to oc¬ 
cur due to differences in self-shielding of the FUV inten¬ 
sity at the dissociation wavelengths for the various CO 
isotopes. The more abundant ^^CO and ^^CO isotopes 
are effectively self-shielded at relatively low extinctions 
(Ay < 1 magnitudes) but because the dissociation wave¬ 
length of C^®0 is slightly shifted compared to ^^CO and 
^^CO, and its abundance is relatively low so it is not 
as effectively self-shielded as the main isotopes, FUV ra¬ 
diation can pen etrate deeper into t he cloud and disso¬ 
ciate the C^®0 (Visser et al. 20091. Since our surveys 
cover relatively high extinctioh regions (Ay > 5 mag) we 
would not necessarily expect selective photodissociation 
to strongly effect the ^^C O to C^^O abundan ce r atios. 
However, ear l ier st udies of Lada et al. (1994) and Shi- 
majiri et al.| (|2014|) suggest that the effects of selective 
photodissociation may be present at much higher extinc¬ 
tions. 


40 



Av (mag) 


Figure 10. The ratio of Wi 3 (^q^ 2 -i) ^Ci® 0 ( 2 - 1 ) plotted 

as a function of Ay- Wi3co(2— i) everywhere > Those 

pixels with Wqi 8 o( 2 -i) ^ plotted as filled circles, those 

originally with Wqi 8 o( 2 -i) < 2c7vy,i8 ^-re set to have W(^i 8 o( 2 -i) 
= 2cr^ 18 to reduce scatter in the diagram and are plotted as trian¬ 
gles. These points are all significant lower limits to the true values 
of the ratio. The white contours show the surface density of pixels 
with detections. At low extinctions (Ay < 15 mag) there is a large 
spread in the observed ratios with the largest values (~ 30-40) 
considerably in excess of the solar abundance ratio of 5.5. At high 
extinctions (Ay > 15 mag) the dispersion in the observed ratios is 
dramat ically less and the ratios are close to the solar value. See 
Section l4.1.5l 


In Figure [TQ| w e plot the ratio of integrated intensities, 
W(^^CO)/w^^®0), as a function of extinction. This 
ratio is proportional to the abundance ratio for the two 
species for optically thin lines. It has been shown to be 
sensitive to the effects of UV radiation on cloud chem¬ 
istry and abundances (e.g. Lada et al.|1994 Shimajiri et 


al.|2014D . For a solar abundance and optical 


y thin lines, 

the ratio should be equal to about 5.5 and one would ex¬ 
pect the points to scatter around this value, independent 
of extinction, provided the cloud is characterized by a 
constant relative abundance throughout and that the ex¬ 
citation of the two species is not too different. Although 
at lower extinctions (Ay <15 magnitudes) there indeed 
appears to be no obvious correlation between this ratio 
and extinction, the observed values scatter over an enor¬ 
mous range (« 0 - 40). The highest values of the ratio 
exceed the solar value by nearly an order of magnitud^ 
In contrast, at high extinctions (Ay > 15 magnitude^ 
the relation is relatively flat with a relatively low dis¬ 
persion, as might be expected for a constant abundance 
ratio, and is centered at a value (approximately 4.5) near 
but somewhat less than the solar value. This behavior 
does not appear to depend on gas temperature (see fig¬ 
ure) and thus is not likely a result of depletion or des¬ 
orption processes on grains. It is more likely indicative 
of an active chemical processing of C^®0 by FUV radia¬ 
tion in regions where Ay < 15 magnitudes. The scatter 
in Figure [T0| is so large at these extinctions because the 
relative abundance of and C^®0 is highly unstable 
possibly due to the stochastic variations in FUV flux re¬ 
sulting, in turn, from such factors as the cloud structure 
and geometry and nature of the external radiation field. 

In Figure [TT] we show a map of the distribution of the 
integrated intensity ratio over the southern portion of 
the CMC. Examination of this map shows systemati¬ 
cally increasing values of the integrated intensity ratio 
at the edges of the cloud with the material nearest the 
star LkHa 101 exhibiting the largest center to edge in¬ 
creases in the ratio. This latter region appears to contain 
a well developed photon dominated region (PDR) while 
the dense cloud material to the north is much more qui- 
escent in this regard. In their study of the Orion cloud 


Shimajiri et al. (2014) reported averaged values of the 
^•^CO to C^®0 integrated intensity ratio ranging from 
about 10 to 17 in the seven PDRs they observed in that 
cloud. Moreover they also found the variations in this 
ratio to be linked to variations in the external FUV ra¬ 
diation field. 

In the inner, higher extinction regions of the CMC, 
where FUV photons cannot penetrate, the abundances 
are more chemically stable and apparently somewhat be¬ 
low the solar value. Our observations do indicate, how¬ 
ever, that FUV photons are present at relatively large 
(projected) depths (Ay ~ 15 magnitudes) i n the cloud 


confirming results from earlier studies (e.g. Lada et al. 


® We note here that in the outer regions of the cloud, C^®0 is 
often undetected while is still quite strong. Because is 

in the denominator, the ratio is very sensitive to noise in the outer 
regions. However, at all locations detections exceed the 5 a 

level. To reduce the scatter in the plot due to non-detections of 
C ^®0 we set the C ^®0 integrated intensity equal to 2 ctvv.is for 
those pixels where W(C^*0) is below 2 cyvis- We represent the 
resulting lower limits on the calculated ratio as triangles in the plot 
(see figure). 
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Figure 11. A map of the to C^®0 integrated intensity 

ratio in the southern portion of the surveyed region. The contours 
show W(^^CO)/W(C^^O) = 5, 10, 15. The largest values of this 
ratio appear in the outer regions where the chemistry is driven by 
FUV radiation that penetrates deeply into the cloud. It is clear 
that this effect is larger in the southern region of the map (i.e., 
S 35°10^ to 35°20^) which is in close proximity to the embedded 
cluster (open white circle) and its UV radiation field than in the 
more remote northern portion of the region. 
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Figure 12. Tex-Ay relation for hot tiles 11-13. The width of 
the bins is 2 mag. The red solid line shows linear fitting on the 
unbinned data points. This shows the correlation between Tex and 
Ay in the region near LkHcn 101. Typical errors in Tex are always 
below 5%. Typical uncertainties in Ay are between 0.05 and 0.1 
magnitudes. See text. 

at high Ay, there is no doubt of a strong trend between 
the two quantities with Tex fx Ay- Interestingly in this 
region the temperature reaches the CO sublimation tem¬ 
perature at roughly Ay = 6-7 magnitudes. Therefore it 
is possible that, even though C^®0 is being evaporated 
from grains at Ay > 6 magnitudes, (where Tex ^ 17 
K) it is also being selectively photodissociated by deeply 
penetrating UV radiation until depths of Ay ~ 15 mag¬ 
nitudes. For larger values of Ay, dust absorption effec¬ 
tively removes the dissociating UV radiation and higher 
gas (dust) temperatures associated with this region lib¬ 
erate more C^®0 from grain surfaces. Indeed, the C^®0 
abundances at these higher extinctions and temperatures 
may reach levels that also allow some C^®0 self-shielding. 
Consequently, C^®0 does not reach the expected levels 
of increased abundance due to evaporation from grains 
until cloud depths are sufficient to shield the molecule 
from UV radiation. Because of its higher abundance, 
^^CO self-shields at much lower cloud depths and the 
effects of its increased abundance, due to evaporation 
from grains, are observed at lower extinctions (and cor¬ 
responding lower temperatures) than C^®0. 


1994| |Shimajiri et al.||2014|). 

We now address the earlier observation that the rela¬ 
tions between the CO abundances and temperature (Fig¬ 
ure]^ show positive correlations but over different tem¬ 
perature ranges. One might have expected both ^^CO 
and C^®0 molecules to show abundances that increased 
with temperature over similar temperature ranges, since 
their binding energies are very similar. However, the dif¬ 
ference may result from the effects of selective UV disso¬ 
ciation of C^®0. 

In the CMC we also find Tex to be correlated with 
Ay in regions where warm gas is present. Figure |12| 
displays the relation between Tex and Ay for the warm 
tiles (11-13) in our map. Although the scatter is large 


4.2. X(CO) on Sub-Parsec Scales 

Even though H 2 is the dominant constituent of molec¬ 
ular clouds, ^^CO emission is the most accessible tracer 
of such gas in the ISM. It is often the only molecular 
tracer that can be readily detected in distant clouds and 
galaxies and consequently has been used to estimate the 
mass of the molecular compone nt of the ISM in these 
systems (e.g. Bolatto et al.|2013). To make such estima¬ 
tions requires knowledge of the so-called X-factor, that 
is, the CO conversion factor, 


Xco = N(H 2 )/W(COj=i_o) 


( 8 ) 


The standard determination of this factor is via use 


of virial theorem techniques (e.g. Solomon et al. 1987). 
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However, the most straightforward way to determine the 
X-factor is by comparing direct measurements of CO and 
Ay on sub-parsec scales in local clouds. Traditionally 
this is accomplished by measuring the slope of the cor¬ 
relation between the two quantities at low extinctions 
(Ay < 5 mag) where the CO appears to be effectively 
thin (e.g. Dickman 1978[ Frerking, Langer, & Wilson 


1982| |Lombardi et al.||2UU6p. We adopt a slightly dif- 
terent approach and use all our observations to compute 
a global or average X-factor for the CMC cloud, that 
is, <Xco(i-o)> = <N(H 2 )>/<W(C 0 )> = 2.53x1020 
cm“2 (K km This value is 25% larger than the 

standard Milky Way value but within the range of values 
deriv ed for individual mo lecular clouds in other studies 
(e.g. Bolatto et al. 2013). Another approach is to de¬ 
rive the on-the-spot Xco, that is, we evaluate Equation 
^ at each pixel in our map and obtain <X(;;;o(i-o)>ots 


0 

= <N(H 2 )/W(C 0 )> = 3.1x102 


(K km s 


The former <Xco> is equivalent to the latter weighted 
by the integrated intensity. 


However, both recent observational studies (Pineda et 
ar]|2008[ 120101 IBieging et al.||2010] [Ripple et al.|| 201 ^ 


al.||2011a|bp 


Lee et alT 2014|) and numerical simulations (IShetty et 
ialbl) have shown that Xco can be strongly de¬ 


pendent on local physical conditions. Since the CMC 
displays a wide range of physical conditions (e.g. tem¬ 
perature, star formation activity and spatially variable 
CO abundances), it would be interesting to measure both 
Xco and its internal variation on sub-cloud, sub-parsec 
scales. 



Figure 13. A pixel-by-pixel comparison between ^^CO(2-l) inte- 
grated intensity and visual extinction Ay for the entire surveyed 
area. Each point represents a 19^' pixel which corresponds to ~ 
0.04 pc at a distance of 450 pc. The red solid line represents the 
standard Milky Way Xco with a slope of 3.3 (K km s~^) mag~^. 
The yellow contours show the surface density of points. Uncertain¬ 
ties are com parable to or less than the size of the symbols. See 
Section EM 


mag ^ (jBohlin et al.||1978 

jRachford et al. 112009 

Pineda 

et ai.|2008p and a constant 

line integrated intensity ratio 

Wco( 2 -i)/Wco(i-o) = 6.7, Xco(i-o) is inversely pro- 


portional to Wi 2 co( 2 -i)/Ay in the plot. In the Milky 
Way disk, the typical value is Xc o(i-o) ~ 2xl02° cm“2 
(K km s“^)“^ (e.g., Bolatto et al.||2013 ). This implies a 


reference value Ay/Wco( 2 -i) ~ 0.30 (K km s” ") " mag. 
We plot a red line indicating this reference value in Fig¬ 
ure [l^ with a slope of Wco( 2 -i)/A.v « 3.3 (K km s“^) 
mag“^. At fixed Ay, a point above this line indicates a 
smaller Xco( 2 -i), and vice versa. 

We can see significant scatter in this plot indicating 
significant spatial variations in Xco( 2 -i) on pixel or sub¬ 
parsec spatial scales. In order to investigate the spatial 
dependence of Xqo( 2 _i), we determined the Wi 2 co( 2 -i)“ 
Ay relation in each i ndividual tile of our CO survey, fol- 
lowiM Section 4.1.2 These relations are shown in Fig¬ 


ure [1^ There are several familiar features of these plots 
that are similar to what we have seen in Figure!^ i.e., 
less scatter in the relations within the individual tiles, 
systematic spatial variations in the Wi 2 ( 2 q( 2 _i) - Ay re¬ 
lation across the cloud, and good correlation between 
high Wi 2 co( 2 -i) and high gas temperature (Tex ^ 15 

K). 

The general similaritjrl in the behavior of Wi 2 co( 2 -i) 
and N(^^CO) with extinction is perhaps surprising given 
that the ^ 2 C 0 emission is very optically thick while 
emission is considerably less thick and mostly optically 
thin. Since N(^^CO) is proportional to Wi 3 co( 2 -i), con¬ 
sider that in this situation Wi 3 co( 2 -i) ~ [J^iTex) — 
Ju{Tbg)] X r^Auia, while Wi 2 co( 2 -i) « [J,^{Tex) - 
JviTbg)] X Aui 2 , where Av is the line width of the corre¬ 
sponding spectrum. Unlike Wi 3 co( 2 -i), Wi 2 (;;;o( 2 -i) has 
no dependence on r or column density, but it is propor¬ 
tional to Tex, the gas excitation temperature. 

Particularly interesting in this regard is the behavior 
of Wi 2 co( 2 -i) in tile 13 near where the embedded cluster 
is heating the cloud. Here, there appears to be a clear 
linear relation between Wi 2 co( 2 -i) and Ay to the high¬ 
est measured extinctions (« 25 magnitudes) despite the 
fact that the ^ 2 C 0 emission in this region of the cloud 
is so saturated that it is self-reversed (see Figure [^. Be¬ 
cause the CO emission is so optically thick, this b^avior 
must be primarily the result of the gradient of Tpx with 


Ay observed in this region of the cloud (see Figure 12) 


This is also evident in Figure [T5| which shows an ex¬ 
panded view of the Wi 2 co( 2 _i)-Xv relation for tile 13. 
We can clearly see a temperature gradient in the distri¬ 
bution of color-coded points from low to high Ay. Both 
W(CO) and T increase by a factor of ^ 5-6 from low to 
high Ay. This indicates that the temperature gradient 
accounts for the bulk of the change in W(CO) with Ay 
in this region. Of course, other factors could contribute 
to the increase in W(CO) such as an increase in velocity 
dispersion with temperature, or even increased [CO] due 
to desorption as with C ^®0 and ^^CO as discussed pre- 


To do this we first performed the traditional pixel-by- 
pixel comparison between Wi 2 go( 2 -i) a nd A y for the 
entire surveyed region, as shown in Figure [Tsl Assuming 
a constant gas-to-dust ratio N(H 2 )/Ay =9.4x102° cm “2 


° By “similarity” we mean the hot tiles 11-14 show ste eper 
slopes than the cold tiles in both Wi2co{2-l) " (Figure |14| 
and N(^^CO) - Ay (Figure!^ relations. The apparent relative dif¬ 
ferences between the the coTa tiles in Figure [^compared to Figure 
|14| is a result of different y-axis scaling. 
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Figure 14. A pixel-by-pixel comparison between integrated intensity and Ay within each tile. All tiles are shown in the same scale 

in both axes. The red solid lines represent the standard Milky Way Xqo with a slope of 3.3 (K km s“^) mag“^. 


viously. We do find some evidence for a marginal (35%) 
increase in the velocity dispersion of the hot gas, but its 
magnitude is far from sufficient to account for the in¬ 
creased (x6) measures of W(CO). A desorption induced 
increase in CO abundance most certainly has occurred 
in the hot gas, and although the ^^CO line is very op¬ 
tically thick, we might expect that additional photons 
could leak out in the line wings and possibly increase 
the measured velocity dispersion. But as we already re¬ 
marked, this effect cannot be significant since the change 


in peak brightness is observed to be comparable to the 
change in W(CO) and substantially larger than changes 
in the velocity dispersion. 

Consequently, in the set of tiles near the massive star 
LkHa 101 (11-14) we conclude that heating of CO en¬ 
hances Wi 2 cof 2 -i) in a manner that induces a positive 
correlation with extinction despite the fact that the emis¬ 
sion is saturated. These observations suggest that heat¬ 
ing by OB stars produces gradients in the excitation tem¬ 
perature of the CO emitting gas such that in warm re- 
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Figure 15. Tile 13 from Figure |14| The black solid line is a 
linear regression. The fitted slope corresponds to X^q = 1-5 x 
10^® cm“^ (K km The red solid line represents the fiducial 

Milky Way value, Xco = 2.0 x 10^*^ cm“^ (K km s“^)“^. 


gions W(^^CO) is proportional to Ay and thus to N(H 2 ) 
over a wide range of extinctions. This, in turn, indi¬ 
cates that a single, meaningful X-factor can be derived 
for these hotter regions using the observations. In Figure 
14l tile 13 shows the best correlation between W(^^CO) 


and Ay. We performed a linear regression for tile 13, 
which is shown in Figure [TSl Our result implies 


-i\-i 


Xqq ~ 1.5 X 10^^ cm ^ (K km s ) 


in this tile over the range 3 mag < Ay < 25 mag. From 
the average of the individual, on-the-spot, Xco values in 
tile 13 we hnd <Xco> = 1-3 x 10^° cm“^ (K km 
in excellent agreement with the value derived from the 
linear regression. This clearly supports the idea of a con¬ 
stant value in the hot gas. The derived value is less than 
both that found for the global average of the CMC and 
the value (2.Ox 10^° cm“^ (K km s“^)“^) usually adopted 
for the Milky Way. 

This result would suggest that in regions heated by OB 
stars (e.g., as HII regions, extragalactic starbursts, etc.), 
reasonably accurate gas masses could be derived using a 
standard X-factor analysis with a single value of the X- 
factor. Unfortunately it is not at all clear how one would 
predict the appropriate value to use in all situations. 

However some clues may be found in the analysis of the 
regions away from OB stars, where there is only ambi¬ 
ent heating of the gas from cosmic rays and the general 
galactic UV radiation held. Here we observe no corre¬ 
lation between Wi 2 co( 2 -i) and Ay at modest to high 
extinctions (i.e., Ay > 3-5 mag). This behavior is due 
to the fact the ^^CO line is almost always saturated at 
high extinctions and the cold gas is very nearly isother¬ 
mal. Consequently, these cold, quiescent regions, are not 
characterized by a single empirical X-factor, instead, the 
X-factor systematically increases with extinction and, in 
the CMC, is almost always larger than the Milky Way 
value. This is shown in Figure |16^ where we plot the 
relation between the on-the-spotX-factor and Ay for 
the regions of the cloud not heated by the cluster (i.e., 
tiles 1-10 and 15-17). As expected there is a clear cor¬ 
relation and a linear least-squares fit to the unbinned 


data hnds that Xco oc N(H 2 )°'^. For the cold tiles we 
find <Xco> = 3.4 x 10^° cm“^ (K km s~^)“^. The 
(constant) X-factor in the hotter regions (e.g., tile 13) 
has a value considerably lower than those in the cold 
regions, suggesting that Xco rnay also increase with de¬ 
creasing gas temperature within the cloud. Such a trend 
is clea rly evident in our observations and is shown in 
Figure 16 d where we plot the relation between the on- 
the-spot Xco and Tex for the entire cloud. A linear 
least-square fit to the (binned) data gives, Xco = 2.0 x 
1020 (t^x/ 10)“°’^ cm“^ (K km s“^)“^. In this instance 
we chose to fit the binned data to give greater weight to 
high temperature measurements. A fit to the unbinned 
data gives a slightly steeper slope (-0.9 vs. -0.7). 

The observed inverse correlation between Xco and Tex 
is potentially interesting. It suggests that differences 
between globally averaged X-factors derived for differ¬ 
ent clouds could be a result of variations in the relative 
amounts of hot and cold gas within the clouds, provided 
their velocity dispersions are similar. Active star form¬ 
ing clouds with OB stars and HII regions, such as Orion, 
might be expected to have average X-factors that are 



cloud complex associated with the HII regions, Xco = 
1.3 X 10^° cm“^ (K km s~^)“^, very similar to what we 
find in the California cloud. For the Taurus cloud, w hich 
has no OB stars or HII regions, Pineda et al.|(2010| find 


Xco = 2.1 X 10^'^ cm ^ (K km s ) ^ from CO ob¬ 
servations of that cloud. A relatively high value (2.54 x 
10^° cm“^ (K km s~ )“^ ) for Xco was also derived from 
Planck observatio ns of high latitude diffuse and pre sum- 
ably cold, CO gas (|Planck Collaboration et al.|2011|). We 
noted earlier that the global X-tactor tor the CIVIC was 
25% higher than the Milky Way value and this is may be 
due to the relatively large tracts of cold regions included 
in the average for the CMC. Of course, variations in ve¬ 
locity dispersions between clouds could also contribute to 
the variation in their observed X-factors. However, we 
would expect similar velocity dispersions to character¬ 
ize gravitationally bound or virialized clouds of similar 
mass and size. Moreover, since the average colu mn den¬ 
sities of molecular clouds ar e also constant (e.g., Larson] 


1981 Lombardi et al. 2010), we would expect <Xco> 
<T>~^ for such ^uds. A s recently pointed out by 


Narayanan & Hopkins (2013), Galactic GMCs are char- 
acterized by very similar physical properties resulting in 
X-factors that span a relatively narrow range of values. 
They suggest that this is a result of the competition be¬ 
tween stellar feedback and gravity that is necessary to 
maintain the relative stability of these objects. Here we 
are suggesting that variation within that narrow range 
of observed Xco may be due to variations in the rela¬ 
tive fractions of hot gas in these clouds. This hot gas, of 
course, is a direct product of the amount of stellar feed- 
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Figure 16. Relations between (a) Xco a-ud N(H 2 ) in cold tiles (excluding tiles 11-14), and (b) Xco and Tex in the entire California 
Cloud. Red (a) and blue (b) lines represent least-square fits to unbinned and binned (every 0.09 dex) data, respectively. The typical 
uncertainties in the X factor are between 5-7%. 


back the clouds are experiencing. Resolved maps of the 
internal X-factor distribution in additional clouds would 
be useful to determine the relative contribution of tem¬ 
perature to the derived values of Xqq within GMCs in 
general and could directly test this idea. 

Absent of such additional data, it is instructive to com¬ 
pare our observations to recent simulations of turbulent 
cloud evolution that resolve synthetic clouds a nd derive 


the internal spatial distribution of the X-factor (Shetty et 
^|2011a|b ). These simulations are thus directly compa- 


rable to observations such as ours. The Shetty et al. sim¬ 
ulations include chemical evolution and radiative transfer 
to enable predictions of CO line emission. However, the 
calculations do not include gas-dust chemical interactions 
and are thus not relevant to our earlier findings or discus¬ 
sion regarding pCO], [C^SO], and W(i3CO)/W(Ci®0). 
But since ^^CO is so optically thick we safely assume 
that effects of depletion/desorption are not important for 
our X-factor determinations and thus comparisons with 
these models are appropriate. Because the Shetty et al. 
simulations do not incorporate stellar feedback, the most 
appropriate comparison will be with our data in the re¬ 
gions of the CMC away from the embedded c luster (all 


tile s exclud ing 11-14). The simulated cloud of Shetty et 


al. (2011b) whose physical conditions match most closely 


Hiose ot typical GMCs is characterized by a global Xco 
lower than we observe in the cold tiles of the CMC map 
(2.2 vs. 3.4 X 10^° cm“^ (K km but closer to 

the Milky Way value. On sub-cloud scales [Shetty et] 
^ (|2011b) find a correlation between Xco and Ay (see 


their figure 4) similar to but somewhat weaker than we 
observe in the CMC. The simulations also predict that 
there should be a correlation between the X-factor and 
the velocity dispersion in the cloud, with Xco cx: <Jyi' 2 - 
In figure [TT] we plot Xco against the velocity dispersion, 
avi 2 in tM cold cloud regiontQ We find no correlation 
between these same quantities contrary to the predictions 


^ o'Vl 2 is the second moment of the ^^CO(2-l) profile = 
( J Tmb(y — Ml) where Ml is the first moment f vTmbdv. o-y 12 

equals the velocity dispersion if line profile is Gaussian. 


of the simulation. Although the simulations of individ- 
ual clouds do not include any systematic temperature 
gradients or increases due to stellar heating, simulations 
were performed to investigate how increasing the over¬ 
all cloud temperature can effect Xco- These simulations 
predict a trend of decreasing X with T, consistent with 
what we observe here. The simulations find Xco (x T“°'^ 
for different temperature clouds which is only somewhat 
weaker than what is observed here, i.e., Xco fx T“°-^ 
within the CMC. Overall, we conclude that the agree¬ 
ment between these simulations and the CMC is mixed 
and it is unclear whether any similarities between the 
simulations and the observations of the California cloud 
are more than coincidental. 
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Figure 17. Xco-ay^i 2 relation for cold tiles, excluding tiles 11- 
14. The yellow contours show the surface density of points. There 
appears to be no strong correlation between these two quantities. 


5. SUMMARY & CONCLUSIONS 


















20 


Kong et al. 


In this study we presented the results of extinction 
and molecular-line mapping toward the most active star 
forming region within the nearby California CMC. The 
relationship between dust and gas properties in this ac¬ 
tive portion of the cloud was investigated with the fol¬ 
lowing results: 

1. The LTE molecular abundances of and C^®0 
in the cloud exhibit significant spatial variations. These 
variations are correlated with the spatial variation in gas 
(and presumably dust) temperature in the cloud. These 
temperature variations are likely caused by the heating 
from the adjacent, partially embedded, LkHa 101 clus¬ 
ter. We argue that the spatial variations in the derived 

and C^®0 gas abundances are largely due to tem¬ 
perature dependent gas depletion/desorption on/off dust 
grains. 

2. The abundance ratio between and C^®0 in¬ 

creases in lower extinction regions of the California cloud, 
suggesting selective photo-dissociation of C^®0 in those 
regions. This is likely caused by ambient UV radiation 
which apparently penetrates relatively deeply (i.e., Ay < 
15 mag) into the cloud, particularly near the embedded 
cluster. There is also evidence that the selective photo¬ 
dissociation of C^®0 can suppress its abundance even 
in warm regions heated by the embedded cluster where 
evaporation off grain surfaces would otherwise lead to an 
increase in the C^®0 abundance at these cloud depths. 

3. Dramatic spatial variations are also observed in 

the relationship between the ^^CO integrated intensity, 
Wi 2 co( 2 -i)) and Ay, particularly at the intermediate 
and high extinctions (Ay > 3-5 mag) primarily probed 
in this study. However, unlike the case for its two rarer 
isotopologues, the variation in ^^CO integrated inten¬ 
sity appears to be a direct result of a spatial gradient of 
the excitation temperature (and gas kinetic temperature) 
with extinction and not the result of any detectable varia¬ 
tion in the abundance due to depletion/desorption 

effects. 

4. We compute the X-factor (=N(H 2 )/W(CO)) for 
each individual pixel in our map and average the results 
to obtain <Xco> = 2.53 x 10^° cm“^ (K km s“^)“^, a 
value somewhat higher than the Milky Way average. 

5. On the sub-parsec scales in the CMC there is no sin¬ 
gle empirical value of the ^^CO X-factor that can charac¬ 
terize the molecular gas in cold (Tk < 15 K) regions with 
Ay > 3 magnitudes. For those regions we find that Xco 
oc Ay^ with <Xco> = 3.4 x 10^° cm“^ (K km s~^)“^. 
We do however find a clear correlation between W(^^CO) 
and Ay in regions containing relatively hot (Tex ^ 25 K 
) molecular gas at Ay > 3 magnitudes suggesting that, 
unlike the cold gas, the warm material may be char¬ 
acterized by a single value of Xco- However in this 
warm gas we find a value for the X-factor, Xco = 1-5 
X 10^° cm“^ (K km s~^)“^, significantly lower than the 
averages for the cold gas, the overall CMC and the Milky 
Way. 

6. Overall we find an (inverse) correlation between 
Xco and Tex with Xco oc Tex~°'^. Such a correlation 
may potentially explain the observed variations in global 
X-factors between GMCs as being due in large part to 
variations in the relative amounts of warm gas heated by 
OB stars within the clouds. 
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